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Recently we have presented a hidden variable model of measurements for a qubit where the hidden 
variable state space dimension is one-half the quantum state manifold dimension. The absence of 
a short memory (Markov) dynamics is the price paid for this dimensional reduction. The conflict 
between having the Markov property and achieving the dimensional reduction was proved in [A. 
Montina, Phys. Rev. A, 77, 022104 (2008)] using an additional hypothesis of trajectory relaxation. 
Here we analyze in more detail this hypothesis introducing the concept of invertible process and 
report a proof that makes clearer the role played by the topology of the hidden variable space. 
This is accomplished by requiring suitable properties of regularity of the conditional probability 
governing the dynamics. In the case of minimal dimension the set of continuous hidden variables 
is identified with an object living an iV-dimensional Hilbert space, whose dynamics is described by 
the Schrodinger equation. A method for generating the economical non-Markovian model for the 
qubit is also presented. 



I. INTRODUCTION 

One of the most peculiar features of quantum me- 
chanics is the exponential growth of resources required 
to define the quantum state of a composite sys- 
tem. It makes the direct simulation of even a handful 
of particles impossible in practice. This growth is due to 
the fact that \ip) contains the full statistical information 
about the probabilities of any possible event, such as the 
joint probability p(si, sjv) of obtaining the outcomes 
{si,...,sat} by measuring the z axis components of N 
1/2-spins. The information in p(s\, sn) grows expo- 
nentially for a given accuracy and considerably exceeds 
the classical information required to specify the actual 
measurement outcome. Since the quantum state is not a 
physically accessible observable, but can be statistically 
reconstructed only by performing many measurements 
on different replicas [l| , it is natural to wonder if this re- 
source excess is strictly necessary to describe the actual 
state of a single realization. The quantum probabilities 
could be reproduced by a hidden variable theory where a 
single system carries less information than the quantum 
state. In such a theory the quantum state \ip) is mapped 
to a probability distribution on a space X of hidden vari- 
able states, i.e., 
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It is clear that the sampling space X can have in principle 
a smaller dimension than the quantum state manifold. 
For example the space of functions on a one-dimensional 
domain is infinite-dimensional and any finite-dimensional 
Hilbert space can be embedded within it. 

In accordance with recent terminology, we will refer to 
the actual state X of a quantum system and the corre- 
sponding space as ontic state and ontological space Q, 
respectively, and name the dimensional reduction of the 
ontological space ontological shrinking. It is interesting 
to note that, in any known short memory (Markov) hid- 
den variable theory, the dimension of the state space is 
never smaller than the quantum state manifold dimen- 



sion. As an example, in the de Broglie-Bohm model the 
wave-function has the role of a field of physical quanti- 
ties and is supplied by additional variables describing the 
particle positions. 

The ontological shrinking has a connection with the 
concept of classical " weak simulation" 0, H| in quantum 
information theory. In a classical "strong simulation" of 
a quantum computer, the goal is to evaluate the measure- 
ment probabilities with high accuracy. This requirement 
is stronger than necessary, since in a real quantum com- 
puter a single run does not give the measurement prob- 
abilities and the output is a precise event. The proba- 
bilities concern the behavior of many experimental real- 
izations. The goal of a classical "weak simulation" is not 
to compute the probability weights, but the outcomes 
in accordance with the weights. There are examples of 
quantum circuits that cannot be efficiently simulated in 
the strong way, but whose weak simulation is nevertheless 
tractable In a hidden variable theory with reduced 
sampling space, the evaluation of the actual dynamics 
of a single realization would require less resources than 
the computation of the quantum state dynamics. Thus, 
the ontological shrinking could offer in a natural way an 
efficient method of "weak simulation" of quantum com- 
puters. 

In the recent years the possibility of a statistical repre- 
sentation of quantum states on a reduced sampling space 
was discussed by various authors [Mil. The problem 
of the smallest dimension of the ontological space was 
posed in Ref. @. It was subsequently proved that the 
ontological dimension cannot be smaller than the quan- 
tum state manifold dimension in the case of a Markov 
hidden variable theory with an additional hypothesis of 
trajectory relaxation Q. We will refer to this result as 
the no-shrinking theorem. Recently we reported an ex- 
ample of hidden variable model of measurements for a 
qubit whose state space is one-dimensional, i.e., smaller 
than the two-dimensional Bloch sphere [l(| • As a conse- 
quence of the dimensional reduction, the dynamics is not 
a Markov process. This counterexample makes evident 
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that the short memory hypothesis is strictly necessary 
for the proof of the theorem in Ref. 7]. In this article, 
we review the one-dimensional model providing a method 
for generating it and analyze in more detail the hypoth- 
esis of trajectory relaxation. We define the concept of 
invertible process and show that it is always possible to 
find a sub-region of a compact ontological space where 
all the processes are invertible. Discarding insignificant 
transient states and considering only invertible processes, 
we present a new version of the no-shrinking theorem 
that makes clearer the role played by the topology of the 
hidden variable space. This is accomplished by requir- 
ing reasonable properties of regularity of the conditional 
probability governing the dynamics and explicitly using 
them in the theorem proof. By the way, it is useful to 
remind that the dimension is a topological property and 
is not defined by the space cardinality. 

In section [TT] we introduce the general framework of a 
hidden variable theory. In section [TTT] the economical on- 
tological model in Ref. is reviewed. The properties of 
regularity and the concept of invertible process are intro- 
duced in section HVl where we also prove the no-shrinking 
theorem and discuss its consequences in terms of resource 
cost. In the same section we show that in any Markov 
hidden variable theory with minimal space dimension 
the set of continuous hidden variables can be identified 
with an object, living in an iV-dimensional Hilbert space, 
whose dynamics is described by the Schrodinger equa- 
tion. In appendix, we report a systematic construction 
method to generate the model in Sec. IIII1 starting from 
a particular form of the probability distribution. 



general also P could depend on the context of the mea- 
surement, thus we introduce additional parameters t in 
the conditional probability, i.e., 



-)• P((p\x,n,r). 



(3) 



The probability of the event \<fi) given \ip) has to be equal 
to the quantum mechanical probability, 



^2 J d M xP(<f)\x,n,T)p(x,n\ip,n) = \ (<j)\ip)\ 



(4) 



Finally, it is assumed that the dynamics at the ontolog- 
ical level is Markovian [l2[ • The ontic state evolves from 
X\ at time t\ to X 2 at time t 2 according to a conditional 
probability K(X 2 , ^l^i^i)- K satisfies the Chapman- 
Kolmogorov equation [12[ and is a delta distribution for 

t 2 = tli 



K{x 2 ,n 2 ,t 1 \xi 1 n u ti) = 5 nu7l2 5(x* 2 - x\). 



(5) 



For a time-homogeneous process the transition probabil- 
ity depends only on the time difference t 2 — t\, 



K (X 2 , t 2 \X u ti) = K(X 2 , t 2 - h\X u 0) 
= K{X 2 \X u t 2 -h) 



(6) 



In order to link the quantum language with the classi- 
cal one, we can label K with the corresponding unitary 
operator U . As for the state preparation and measure- 
ment, in general the conditional probability can depend 
on additional parameters Xi the transformation context, 



II. GENERAL FRAMEWORK 



U^KfijXilXi) 



(7) 



In a general hidden variable theory the quantum state 
is translated into a classical language by replacing it with 
a probability distribution on a sampling space X of ontic 
states. We assume that the ontological space is a M- 
dimensional manifold described by an M-tuple x of con- 
tinuous variables and a possible discrete index n. The 
mapping ([1]) is not the most general, since the proba- 
bility distribution could depend on the preparation con- 
text. For example, a 1/2-spin can be prepared in the 
up state by merely selecting the beam outgoing from a 
Stern-Gerlach apparatus or can be prepared in the state 
1 1) + I 1} an d then suitably rotated. In order to account 
for this possible dependence, we add suitable parameters 
?7 that identify the preparation context llj, i.e., 



(2) 



The set of ontological variables X — (x, n) contains the 
whole information about a single realization, thus the 
probability of any event is conditioned only by it. In 
particular, for the measurement of the trace-one projec- 
tor |(/>)(<^|, it is assumed that there exists a conditional 
probability P for the event \(j>) given the ontic state X. In 



Indeed an operator U can be physically implemented in 
different ways. For example, the spin rotation e _itcrx can 
be performed directly rotating along the x axis or imple- 
menting the three-step rotation e~ l T ay e~ lt(7z e*^ * along 
the z and y axes. These two schemes are physically dif- 
ferent and not necessarily described by the same condi- 
tional probability. By the way, it is useful to note that 
Kjy x (X\X) for U — 1 is not necessarily the delta distri- 
bution S(X — X) for all the contexts. For example the 
identity evolution can correspond to physically perform- 
ing the three step rotation e ~ lJ i cry e~' lt ' 7 'e l ^ cry , followed 
by the inverse transformation e lt<Jx . The overall opera- 
tion is not equal to doing nothing and does not necessar- 
ily correspond to a delta peaked conditional probability. 

If the quantum state evolves to U\ij)) = the 
associated probability p(X\tp,r]) evolves to 



dYK jX\Y)p(Y\^,r,). 



(8) 



Some regularity properties of Kjy ^ will be discussed in 
Sec. IIVI Any short memory hidden variable theory has 
this general structure. 
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III. ECONOMICAL ONTOLOGICAL MODEL 
OF MEASUREMENTS 

In this section we show that if the Markov condition is 
not assumed, then the ontological space can be smaller 
than the quantum state manifold. This goal is achieved 
by explicitly providing a one-dimensional hidden variable 
model for a qubit [10(. Its systematic construction is 
discussed in the Appendix. 

For the moment we introduce a model working only 
for a subset of preparation states. The extension to the 
whole quantum state manifold will be discussed later on. 
The ontological space is given by a continuous real vari- 
able x and a discrete index n that takes the two val- 
ues and 1. It is convenient to represent the quantum 
state \ip) and the event \<f>) by means of the Bloch vectors 
v = (i/j\a\ip) and w = (<f>\a\<f>), where a — (a x ,a y ,a z ), &i 
being the Pauli matrices. 

The probability distribution associated with the state 
v is 



p(x, n\v) 



5 n<1 S(x-6), (9) 



where 9 and tp are respectively the zenith and azimuth 
angles in the spherical coordinate system 



v x 

Vy 



sin v cos ip, 
sin 9 sin tp, 
cos 9. 



(10) 



Thus, when the quantum state v is prepared, the index n 
takes the value or 1 with probability sin 9 or 1 —sin 9 and 
the continuous variable is equal to the zenith or azimuth 
angle according to the value of n (see Fig. [TJ. 



n=0 
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FIG. 1: One-dimensional ontological space {x,n} at left. Ac- 
cording to Eq. @, each point of the Bloch sphere is asso- 
ciated with a probability distribution on {x, n} having two 
delta peaks. 

The conditional probability P(w\x,n) for an event w 
with w, > is defined as follows 



p(w\ x , o) = i + w * cosx + w y sin x - vlEM, (ii) 



P(w\x,l) 



wl sin x + w z cos x — 1 



2(1 — sin x) 



(12) 



The events w with w z < correspond simply to the 
non-occurrence of the events —w with w z > 0, i.e, 
P(— w\x, n) = 1 — P(w\x, n). 

It is easy to prove that these probability functions fulfil 
the condition (J4}, that is, 

P(w\tp, 0) sin 9+ P(w\9, L)(l-sin0) = i (1 + w ■ v) (13) 

As shown in Ref. [10], the conditional probabilities 
P(w\x, 0) and P(w\x, 1) are always smaller than or 
equal to 1, but P(w\x, 1) is positive only if 9 < 9q = 
arccos(|) ~ 53.13 degrees. 

Thus, the model in this form works only for a set of 
prepared states whose Bloch vector lies inside a cone with 
aperture 29q, the z axis being the symmetry axis. Since 
the positive region has non-zero measure in the quan- 
tum state manifold, it is possible to extend the model to 
the whole Bloch sphere by covering the manifold with a 
sufficiently large number of patch regions with different 
coordinate systems and enriching the ontic state with a 
finite quantity of information. This was accomplished in 
Ref. [10( by adding a discrete index m taking 12 possible 
values, labelling 12 different regions of the Bloch sphere. 
See the referred paper for further details. 

This is a concrete example of ontological shrinking, 
where the ontic state space is smaller than the quantum 
state manifold, which is in this case the two-dimensional 
Bloch sphere. It is a remarkable fact that a single real- 
ization of the ontic state {x, n} contains less information 
than the quantum state. The whole information on \ip) 
is contained in the probability distribution p. 

The absence of a short memory description of the dy- 
namics is the price paid for the dimensional reduction. 
Indeed it is impossible to associate a positive conditional 
probability with each unitary evolution, that is, it is im- 
possible to satisfy the identity 

p(x,n\Utp) = 2_. / dxKjj(x,n\x,n)p(x,n\-ip) (14) 



with the probability distribution in Eq. ([5]), apart from 
the unitary evolution e IttT3 . Indeed in the other cases the 
dynamical equation of p(x,n\ip) is in general nonlinear. 
Let us consider the unitary evolution with the Pauli ma- 
trix a y as generator. The dynamical equations of the 
Bloch vector are 



fr =° 

ov z _ 

at ~ 



Vz, 



(15) 



which correspond in spherical coordinates to 

|# = — cot 9 sin ip 
§ =cos</?. 



(16) 



Let K(x,n\x,n,i) be the transition probability associ- 
ated with the evolution in Eq. (fT5")) . we have from Eq. (Ti"4")) 
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that 



dp(x, n, t) 
di 



E 



.OK . . . . 

ax-g—(x, n\x, n, t)p(x, n), 



(17) 



that becomes by means of Eq. ((9|) 



sin^o^f^iff + (l-sin< 



r s(x-e) ae I 
J ™. 1 at" 1 " 

cosfl [<y„, «(x - v) - -°)]w = 

^-(x,n\tp,0,t) sm8 + ^(x,n\9,l,t)(l - sin ff 

In particular, for n — 0, and using Eq. (1161) . 



(18) 



cos 



dip 



sin <£> + cos 06 (x — tp) cos tp 



^f(x,O\ip,O,t)sm0 + ^f(x,O\e,l,t)(l-smt 



(19) 



Dividing both sides by sm9 and differentiating with re- 
spect to 8, we obtain that 



cot g2il^—!fl shop + cot 6*<5(x — tp) cos </? 



(20) 



This equation is not satisfied by any function K, since 
the left-hand side is a function of both 9 and tp, whereas 
the right-hand side depends only on 9. Thus, the dy- 
namical equation of the probability distribution @ is 
non-Markovian. 



IV. ONTOLOGICAL SHRINKING AND 
MARKOV PROCESSES 

In this section we will prove that the ontological shrink- 
ing is in contradiction with a Markov dynamics. This 
purpose is achieved by means of a very reasonable hy- 
pothesis concerning the support of the conditional prob- 
ability K (X2I-X1, t). It will be introduced in the following 
subsection. In subsection IIV Bl we will introduce the con- 
cept of invertible process and will show that it is always 
possible to find a sub-region of a compact ontological 
space where all the processes are invertible. Then we 
prove the no-shrinking theorem in subsection IIV CI In 
subsection IIVD1 it is shown that in the case of minimal 
ontological dimension it is possible to identify the set 
of continuous ontological variables with a vector of the 
Hilbert space whose dynamics is given by the Schrddinger 
equation. In the last subsection we discuss the conse- 
quences of the no-shrinking theorem in terms of resource 
cost. 



A. Conditional probabilities and associated unitary 
operators 

For the sake of simplicity, from now on we will omit 
without loss of generality the discrete index n in the def- 
inition of the ontic state and assume that the ontological 
space is a difierentiable manifold whose points are identi- 
fied by the vector x. When the local Euclidean structure 



is not required, we will use the more generic symbol X 
to indicate the ontic state. 

In order that a transition probability describing a 
Markov process makes physically sense, it has to satisfy 
some conditions of regularity. For our purpose, it is suf- 
ficient to require a very weak condition. As a reasonable 
hypothesis, we assume that there exists a subset of the 
support of the conditional probability K[X2\Xi,t) that 
changes smoothly with respect to the condition X\ and 
the evolution time t. Let us refine this statement in a 
more precise way. 

Property 1: given a time-homogeneous process X\ — > 
Xi with non-zero probability K(x2\x\,t), there exist an 
M x M matrix A and a vector a such that K(x2 + \6x"i + 
a6t\xi + 6x1, t + 6t) =/= for any infinitesimal variation 
6x\ and 5t. The matrix A and the vector d? are functions 
of t and the process Xi — > Xi. 

It is important to note that this property is fulfilled 
for the very large class of Markov processes that involve 
drift, diffusion and jumps (they are discussed for example 
in Ref. J.12|). Furthermore, the matrix A and the vector a 
are not necessarily unique. For example in the case of a 
stochastic process any choice of A and a is suitable, since 
the conditional probability is a multidimensional smooth 
function whose support is the whole manifold. For a 
pure deterministic process both A and a are unique. The 
latter gives the drift velocity of the ontic state. In the 
case of pure jumps with finite transition probability, one 
choice is a — and is unique if the jump distance cannot 
be arbitrarily small. One can imagine more complicated 
cases that involve diffusion in some direction, drift in the 
other ones and jumps, however also in these situations 
the stated property is fulfilled for some A and a. 

Property 1 strictly depends on the topology of the on- 
tological space since it involves its local Euclidean struc- 
ture. Using this structure in the proof of the no-shrinking 
theorem is fundamental because the space dimension is 
not a property defined merely by the cardinality of the 
space. 

The quantum unitary evolution is a continuous pro- 
cess and corresponds to a trajectory in a Lie group man- 
ifold with the Hamiltonian as generator. For a time- 
homogeneous process, the unitary operator of the evolu- 
tion has the form e~ ztH , where H is the transformation 
generator and t is the evolution time. In general, it is 
not possible to directly implement every generator of the 
su(N) algebra by means of a purely time- homogeneous 
process. In practice, only a small set of evolutions (build- 
ing blocks) can be directly generated by a physical pro- 
cess. The other evolutions are obtained by suitably con- 
catenating the building blocks. Given any generator G, it 
will be possible in principle to experimentally implement 
a unitary evolution U\, an physically attainable time- 
homogeneous process e~ ltH and another unitary evolu- 
tion U2 such that 



U 2 e 



-itH 



-itG 



(21) 
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The operators C/j can be constructed by suitably con- 
catenating physically attainable processes. This allows 
us to associate with any unitary operator e~ ttG a con- 
ditional probability K(X 2 \Xi,t) satisfying Property 1. 
Note that t in K does not correspond in general to the 
evolution time of a purely time-homogeneous process. It 
is the evolution time of only a part of the overall pro- 
cess involving also the transformations Ui. This implies 
that for t — the conditional probability K(X 2 \Xi,t) is 
not necessarily a delta distribution. Furthermore, we can 
add a shift time to to t and absorb the extra-term e~ lt " H 
in Eq. (f2~Tj) into the operators U\ and U 2 . Thus, it is 
neither necessary to require that i is a positive quantity, 
as in the case of stochastic Markov processes where the 
propagation kernels are elements of a semi-group. For 
our purpose, it sufficient that K(X 2 \X\,£) is defined in 
a neighborhood of t = 0. 

Let {Gi} with i = 1, ...,D be a set of D = N 2 - 
1 generators of the Lie algebra. The unitary operators 
e -itGi are assoc i a ted with the conditional probabilities 
Ki{X 2 \X x ,t). 

Any unitary evolution U in a region around the iden- 
tity can be constructed in the following way, 



U(ti, ..,*£>) = e 



_ c -it 1 G 1 



^ — it D Go 



D 

n« 



-itiGi 



(22) 



where the variables t< parametrize the SU (N) manifold. 
Let the linear operator 



P{X) 



dYKi(X\Y, t)p(Y) 



be denoted by Ki(t). The overall evolution U(t) is asso- 
ciated with the conditional probability 



K(X 2 \X 1 ,t 1 ,-,t D ) 



{X 2 \X X ), (23) 



the product order is such that the sum index grows from 
left to right. 

It is easy to prove a general property of the conditional 
probability K(X 2 \Xi , t) . 

Property 2: suppose that K(x 2 \xi,t) is different from 
zero for some process X\ — > X 2 , then there exists aMxD 
matrix f] such that, for any small variation of the time 
parameters 



St = 



( Sh 

St 2 

\St D 



the conditional probability K(x 2 + fjSi\xi,t + St) is differ- 
ent from zero. In the following, fj will be called the shift 
matrix of K(X 2 \Xi,t). As with the matrix A and the 
vector a, it is a function of t and the process Xi X 2 . 



Proof: first, we consider the concatenation K\ 2 of two 
conditional probabilities, that is, 

K 12 {X 2 \X u t u t2) = J dZK x {X 2 \Z,t x )K 2 {Z\X x ,t 2 ). 

(24) 

Let X\ — > X 2 be a process with non-zero transition prob- 
ability Ki 2 (X 2 \Xi,ti,t 2 ). There exists a state Z such 
that K\{X 2 \Z,ti) and K 2 (Z\Xi,t 2 ) are different from 
zero. Property 1 implies that there exist two vectors <Ji 
and o?2 and a matrix Ai such that 

K x (y + \\d 2 5t 2 + 6tx6ti\z + a 2 St 2l ti + Sti) ^ 0, 



^2(^+32^21^*2 + ^2) + 0. 

Thus, the conditional probability K\ 2 (x 2 + fj5t\xi^ti + 
St±, t 2 + St 2 ) is different from zero, where the columns of 
the M x 2 matrix fj are c?i and Aid^- The property can be 
proved by induction for any concatenation of processes 
K t . □ 

There is an important direct consequence of this prop- 
erty. 

Lemma 1: A process X± — > X 2 with non-zero proba- 
bility K(X 2 \Xi,Q) is associated with a Z? s -dimensional 
manifold of unitary transformations U(St), where D s > 
D — M, M being the ontological space dimension. The 
manifold is identified by the M equations 



D 



= 0, with i = 1, 

3=1 



..,M, 



(25) 



where fj is the shift matrix of K(X 2 \X\, 0). In particular, 
D s = D—Nj, Nj(< M) being the number of independent 
equations in the constraints (|25[) . 

Proof: because of Property 2, there exists a M x D 
matrix fj such that K(x 2 + f)5i\xi,6i) is different from 
zero, i.e, x\ — > x 2 + fjdt is a process associated with the 
unitary evolutions U(St). In particular the submanifold 
of unitary evolutions with fjdt = is associated with the 
process X\ — > X 2 . Its dimension D s is equal to D — Ni, 
Ni being the number of independent equations in fjdt = 
0. Since A/ < M, we have D s > D — M . 

The actual value of D s depends on the number Nj of 
independent constraints in the vectorial equation fjSt = 
0. For example, D s is equal to N 2 — 1 for fj = 0, which 
is the case if Ki(y\x,t) are M-dimensional diffusive pro- 
cesses. 



B. Set S(X) and its symmetry property 

In an ontological model the quantum state is associ- 
ated with a probability distribution p(X\ip,r]) according 
to the mapping © . It is useful to introduce the following 
definition of the set S(X). 

Definition 1: a vector \ip) of the Hilbert space is in 
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S(X) if and only if there exists a context n such that the 
probability p(X\t/j,rj) is different from zero for the state 
X. 

In other words, the set S(X) contains every quantum 
state that is compatible with the occurrence of the ontic 
state X. As discussed in Ref. Q, the set <S cannot lose 

vectors along its evolution. More precisely, if X -H» Y is a 
non-zero probability process associated with the unitary 
evolution U, then US(X) C S(Y). This is a direct con- 
sequence of the definition. Indeed, if \ip) is in S(X), then 
there exists a probability distribution associated with \ip) 
such that X is in its support. Since X — > Y is a non-zero 
probability process, then Y is in the support of a prob- 
ability distribution associated with the evolved quantum 
state U\tp), that is, \ip) E S{X) U\ip) E S(Y). The 
opposite implication E S(X) •<= U\ip) E S(Y) is not 
trivially satisfied. Thus, the set S cannot lose vectors, 
but in principle it could grow acquiring vectors. 

The opposite implication can be deduced by assuming 
that each process is invertible. 

Definition 2: a non-zero probability process X\ — ^ X2 
is said to be invertible if there exists a unitary operator 

U 2 such that X 2 ^4 X 1 is a non-zero probability process. 
Property 3: every process is invertible. 

Note that the operator U 2 associated with the inverse 
process is not required to be necessarily the inverse of 
U\ . Property 3 is very reasonable and is satisfied by any 
known hidden variable theory. Indeed a state connected 
to other states by means of a non-invertible process would 
be only transient and could be safely eliminated by the 
theory. The fact that the transient states are insignificant 
is made clearer if it is assumed that the ontological space 
is compact. 

Proposition: if the ontological space is compact, then 

the processes in any series X\ ^4 X 2 -4 X3 ^4 ... become 
closer and closer to being invertible. 

Proof: a metric ontological space is compact if it is 
closed and bounded. Let Br(X) be M-dimensional balls 
with radius R and center X. Suppose that for any R 

there exists a series X\ -4 X 2 — I X3 -4 .. . with an infinite 
number of processes that take a state X away from its 
ball Br(X). Then there exists a subseries Y\ — y Y 2 — > 
Y3 — > ... where Y m £ Bji(Y n ) for every m > n, but 
this is impossible because the space is bounded. Along 
the series, the processes connecting an element to the 
following ones become in fact closer and closer to being 
invertible. □ 

Discarding insignificant transient states and taking for 
granted Property 3, we can prove the second lemma. 

Lemma 2: Assuming Property 3, if X\ -4 X 2 is a non- 
zero probability process, then U\S{Xx) = S(X 2 ). 

Proof: the process X\ ^4 X 2 is allowed, thus 

U 1 S(X 1 )CS(X 2 ). (26) 

Since there exists an operator U 2 such that X 2 -4 X\ is 



a non-zero probability process, then we have also that 
U 2 S{X 2 )CS(X 1 ). (27) 

These two relations imply that 

US(X 2 ) C S(X 2 ), (28) 

where U = U\U 2 . By iteration we obtain that 

U n S(X 2 ) C S(X 2 ), (29) 

for any integer n. For a finite dimensional Hilbert space, 
it is always possible to find an integer n such that 
U n is very close to the inverse operator U^ 1 . Thus, 
U- 1 S(X 2 ) C S(X 2 ), that is, 

US{X 2 ) D S(X 2 ). (30) 

From inclusions (|28I30[) we have that 

US(X 2 ) = S(X 2 ). (31) 

Applying the operator U\ to both sides of inclusion (|27|) 
and using Eq pip , the inclusion 

S(X 2 ) C UiSiXt) (32) 
is obtained. This relation and inclusion (|26[) imply that 

f/ 1 5(X 1 )=5(X 2 ) (33) 
and the lemma is proved. □ 

Note that, for symmetry reasons, a similar equation holds 
also for the inverse process, that is, 

U 2 S{X 2 )=S{X 1 ). (34) 

In Ref. 0] we proved the following property for S(X). 
Lemma 3: the set S(X) cannot contain every vector 
of the Hilbert space. Equivalently, the set S(X) is not 
invariant with respect to the group SU(N). 

Proof by contradiction: suppose that S(X) contains 
every vector of the Hilbert space, then it contains in 
particular also two orthogonal vectors. This means that 
there exist two overlapping distributions associated with 
two orthogonal quantum states. But this is impossible 
because two orthogonal states can be perfectly discrimi- 
nated by a measurement 0, [ill ] . Indeed the probability 
of obtaining \ip) given \tp) is 

J dXP^\X,T)p(X\^, v ) = |(V#)| 2 = 1. (35) 

This implies that P(tj)\X, r) is equal to 1 in the support 
of p(X\ip,r]). However the probability of \ip) given an 
orthogonal state \if>±) is 

J dXP(il>\X,T)p(X\il> x ,fj) = \(^±)\ 2 =0. (36) 
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This implies that p(X\ip±,r]) cannot be different from 
zero if X is in the support of p{X\ijj,rj), where 
P(ip\X,r) = 1. □ 

Finally, we enunciate the last lemma. 
Lemma 4: Let G be a Lie subgroup of the group SU (N) 
acting on a iV-dimensional Hilbert space T~L and S be a 
set of vectors in H. If the manifold dimension of G is 
larger than (N — l) 2 and S is invariant with respect to 
G, then S contains every vector of the Hilbert space, that 
is, S is invariant with respect to SU(N). 

Proof: any compact Lie group and their linear repre- 
sentations on C N are well-known. One can check that, for 
N ^ 4, the proper Lie subgroup of SU(N) with largest 
manifold dimension is SU(N — 1) x £7(1). Its dimension 
is equal to (N — l) 2 . Thus, if the dimension of G is larger 
than (N — l) 2 , then G is in fact SU(N) and S contains 
every vector. In the special case N = 4, the symplectic 
group Sp(2) is the subgroup of 5i7(4) with largest di- 
mension. A set of generators in the representation space 
<C 4 is 



with i = 1,2,3, a i being two sets of Pauli matrices 
acting on the tensor space C 2 <g> C 2 = C 4 . It is easy 
to check that these generators form the basis of a Lie 
algebra with dimension equal to 10, which is larger than 
(N — l) 2 = 9. However also in this case if the set S is 
invariant with respect to Sp{2), then it contains every 
vector of the Hilbert space. Indeed, it is possible to show 
that the only orbit of Sp(2) is the whole Hilbert space, 
the orbit of a subgroup being the set of states connected 
by means of some subgroup element. It is sufficient to 
prove that any vector is connected to | f) | f) . A generic 
vector \ip) has the form 

|V) = cos 1 t)|0i > +sin0| |)|02 >, (38) 

where |0i) and \<f>2) are two-dimensional vectors. This 

vector is connected through the generators to a vec- 
tor with the form 

= cos6\ t)| t> +sin0| |)|0 2 > . (39) 

Through the unitary operator e 2 +92<7 2 ) anc j 

a suitable choice of 6^, it is possible to connect \ijj) to 

|V) =cosfl| t)| t> +e iv sm9\ |)| |> . (40) 

This last vector is connect to | |)| t) through the gen- 
erators and {^ ] (t {2) + a 2 1) & {2) )/2 = -ia { + ] a {2) + 

ia_a_, a± being the raising/lowering operators. □ 

This theorem implies that the minimal number of pa- 
rameters required to define the orientation of a set S can- 
not be smaller than 2N — 2, apart from the trivial case 
of a completely symmetric set S, which does not have 
an orientation. This property was intuitively introduced 
in Ref. 0, where we noted that a set with the highest 
symmetry SU(N — 1) X U(l) requires 2N — 2 variables 
to specify the orientation of its symmetry axis. 



C. No-shrinking theorem 

At this point we have sufficient tools to prove the no- 
shrinking theorem. We will show that if the ontologi- 
cal space dimension would be smaller than the quantum 
state manifold dimension, then the set S(X) would be 
invariant with respect to the group SU(N), in contradic- 
tion with lemma 3. 

Theorem: If the dynamics in a Markov ontological 
theory satisfies Property 1 and all the processes are in- 
vertible, then the ontological space dimension M is not 
smaller than the quantum state manifold dimension, that 
is, M > 2N — 2, N being the Hilbert space dimension. 

Proof by contradiction: suppose that M < 2N ~2. Let 
us consider the unitary evolution U(t) defined in Eq. (|2"2")l 
and the associated conditional probability K given by 
Eq. (|23p . We have from lemma 1 that a process X — > Y 
with non-zero probability K(Y\X,0) is associated with 
a D s -dimensional manifold of unitary evolutions, say U, 
with 

D s > N 2 - 1 - M > N 2 - 1 - 2N + 2 = (N - l) 2 . (41) 

The manifold contains the identity. As a consequence, 
the set S(X) evolves, for any unitary evolution in U, to 
the same S(Y). By means of lemma 2 we have that 

US{X) = S(Y), VU e U. (42) 

Since leU, 

S(X) = S(Y). (43) 

Thus, 

US(X)=S(X), VC/GW, (44) 

that is, S(X) is symmetric with respect to a group whose 
generators are algebraically generated by D s generators. 
The manifold dimension of this group is equal to or larger 
than D s > (N — l) 2 . Because of lemma 4, S(X) is invari- 
ant with respect to SU(N), but this is in contradiction 
with lemma 3. We conclude that the ontological space 
dimension M cannot be smaller than the quantum state 
manifold dimension 2N — 2. □ 



D. Minimal dimension and SchrSdinger equation 

It is possible to show that in the case of minimal on- 
tological dimension the set of continuous variables can 
be identified with a vector living in the Hilbert space, 
whose dynamics is described by the Schrodinger equa- 
tion. It does not necessarily coincide with the quantum 
state, but moves rigidly with respect to it. 

Suppose that the ontological dimension is minimal, 
that is, M = 2(N — 1). From lemmas 1-4 we have 
that the set S(X) has to be symmetric with respect to 
a subgroup of SU(N) with manifold dimension equal to 
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N 2 -l-M = {N-l) 2 . The subgroup is SU(N-l)xU(l). 
Thus, the set has a symmetry axis \4>), which is a vector 
in the Hilbert space and identifies the orientation of S. 
Since each ontic state X is associated with a set S(X), 
we have the mapping 



(45) 



Because of Lemma 2, in a process associated with the 
unitary evolution U(t), the set S rotates rigidly to U(t)S. 
In particular, its symmetry axis satisfies the Schrodinger 
equation 



(46) 



where H(t) is the Hamiltonian that generates U(t). 

The mapping (|45[) is surjective, but is not necessarily 
injective. However, since the number of continuous on- 
tological variables is just sufficient to label \4>) up to a 
global phase, an additional discrete index is sufficient to 
make the mapping bijective, 



(47) 



Thus, the ontic state in a minimal Markov theory is 
identified with a vector in the Hilbert space and a possi- 
ble additional discrete index. The dynamics of the vector 
is given by the Schrodinger equation. \4>) is equal to the 
quantum state if the set S(X) contains only one state 
that coincides with the symmetry axis. This is the case 
in a wave-pilot theory such as the de Broglie-Bohm me- 
chanics. 

It is interesting to note that requirements (|4^1 |4T|) arc 
fulfilled by the Kochen-Specker model [13|, where the 
ontic state is identified by a Bloch vector having the same 
dynamical equation of the quantum state. 



E. Resource cost with round off error 

We have proved that in any hidden variable theory 
with short memory dynamics the dimension of the state 
space cannot be smaller than the quantum state manifold 
dimension 2N — 2. There is a relation between dimen- 
sion and resource cost required to identify an ontic state. 
By resource cost we mean the quantity of information 
required to identify an ontological state. Obviously the 
information carried by an ontic state is infinite, since the 
ontological space is continuous. Thus, the resource cost 
as a function of the dimension makes sense only in pres- 
ence of a fixed round off error of the continuous variables. 

One can find a scaling law between dimension and re- 
source cost in the following way. Suppose that the ontic 
state is identified by M continuous variables Xi, which 
arc defined in a finite interval. With a suitable rescaling, 
we can assume that they run in the interval between 
and 1. The probability P{£ \x) of an event E is condi- 
tioned by the ontic state x. It is assumed that P{£ \x) is 



a smooth function of x, but the discussion could be ex- 
tended to the case of a finite number of discontinuities. 
Let gi be the mean magnitude of \d Xi P(£\x)\. 

Let each continuous variable be discretized by rii 
points. This introduces a round off error AE in P that 

) 2 for sufficiently large values of 



Y,tLi(9*/ n * 



scales as 

rii. The number of bits required to identify the ontic 
state on the lattice is proportional to the information 



A I 



i = $> g 



(48) 



For a fixed I, using the Lagrange multiplier method we 
find that the optimal choice of rii that gives the smallest 
error is 
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(49) 



where g = (\\i9i) 1 ^ M * s tne log-average of gi. Thus, we 
have the scaling law 



AE ~ gM*e 



(50) 



The error exponentially decreases with X at a rate in- 
versely proportional to M. For a fixed error, the infor- 
mation scales as 



2 ~ M log 



Mig 
AE 



(51) 



The no-shrinking theorem states that M is not smaller 
than the quantum state manifold dimension. Thus, in a 
composite system, I grows at least exponentially in the 
number of parts (for example the number of qubits in a 
quantum computer). 

It is possible to estimate a lower bound for g. Since the 
probability of an event goes from to 1 and xt rambles 
about the interval [0:1], it is reasonable to assume that 
gi cannot be smaller than a value around 1, that is, 
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> 1. 



Thus the linear growth of I with respect M cannot be 
mitigated by an exponential decrease of g. 



V. CONCLUSIONS 

We have shown that in any hidden variable theory with 
a short memory dynamics the ontological space dimen- 
sion cannot be smaller than the quantum state mani- 
fold dimension. Thus, like the quantum state, the ontic 
state necessarily carries for a given accuracy an amount 
of information that grows exponentially in the number of 
subsystems. In comparison with Ref. Q, we have pro- 
vided a better justification of the relaxation hypothesis, 
showing that in a compact ontological space it is possible 
to find a trajectory along which the system unavoidably 
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converges towards a region where all the processes are 
invertible. Furthermore we have presented a hidden vari- 
able model of measurement for a qubit whose ontological 
space is one-dimensional, which is one-half the dimension 
of the Bloch sphere. The corresponding dynamics is not 
Markovian, in accordance with the no-shrinking theorem. 

This model provides a counterexample making evident 
that the hypotheses of short memory is strictly neces- 
sary to prove the no-shrinking theorem. By dropping it, 
we have shown that a single realization can carry less 
information than the quantum state. More drastically, 
we could drop the causality hypothesis. Indeed the non- 
causality is implied also by the Bell theorem and the 
Lorentz invariance of the hidden variable theory. Thus, 
there are two signs that point to the same direction, that 
is, the rejection of causality at the ontological level. The 
possibility of an ontological shrinking for a general N- 
dimensional Hilbert space in a theory without a Markov 
dynamics is an open question whose answer could provide 
a deeper understanding of the computational complexity 
in quantum mechanics. 

I wish to thank J. Wallman and R. W. Spekkens for the 
careful reading of the manuscript and useful suggestions. 
Research at Perimeter Institute for Theoretical Physics is 
supported in part by the Government of Canada through 
NSERC and by the Province of Ontario through MRI. 



Appendix A: Deriving the economical model for a 
qubit 

In this appendix we will present a systematic method 
to generate the model reported in Sec. [TTJ We will con- 
sider the class of probability distributions with the form 



p(x, n\v) = r(n\v)8 [x - f n {v)} 



(Al) 



where the pair (x e R,n G {0, 1}) is the ontic state and 
v = (ip\a\i/j) is the Bloch vector corresponding to the 
quantum state \ip). The quantities f n (v) are real func- 
tions. The normalization of the distribution gives 



r(0\v) + r(l\v) = 1. 



(A2) 



When a system is prepared in v, at the ontological level 
there is a probability r(0\v) [r(l\v)] that the discrete in- 
dex takes the value n = (n = 1). Correspondingly, the 
continuous variable takes with certainty the value fo(v) 
[fi (v)]. This class is the simplest one that fulfils some 
requirements. As noted in Subsec. HVBl two orthogo- 
nal states cannot have overlapping probability distribu- 
tions, implying that the support of p(x, n\v) cannot be 
the whole ontological space. This rules out the class of 
smooth analytical distributions and leads us to consider 
the probability distributions with zero-measure support 
as simplest case. In particular, the distributions with a 
two-point support are parametrized with 3 real param- 
eters [the positions f n (v) of the points and the relative 
probability weight r(0|u) — r(l|u)], which are sufficient in 
order to cover the two-dimensional Bloch sphere. 



For our purposes, it is convenient to define the vari- 
ables x n = f n (v) and to use them to parametrize the 
quantum state. Thus, the probability distribution be- 
comes 



p{x, n\x , x\) = r(n\x , xi)S (x - x n ) . 



(A3) 



Let \<j)) be the event of some projective measurement. As 
for the preparation state, we introduce the Bloch vector 
w = (4>\a\<f)) to label the event \cf>). P n (w\x) is the con- 
ditional probability of the event w given the ontic state 
(x,n). Equation ([4]) becomes 



l 

E 

n=0 



dxP n (w\x)p(x, n\x) 



1 + w ■ v(x) _ 



S, (A4) 



where v(x) gives the Bloch vector v as a function of the 
parameters {xq, x\) = x. The quantity l/2(l+w-v) is the 
Born probability |(0|-0)| 2 in terms of the Bloch vectors. 

Using Eq. (|A3|) . Equation (|A4[) becomes after integra- 
tion 

P o (w\x o )r(0\x o ,xi) + Pi(w\xi)r(l\x ,xi) = S. (A5) 

Note that the conditional probability P n (w\x n ) in this 
equation depends only on the variable x n . Thus, 
in order to find it, we could consider the four pairs 
(x ,xi), (afo,yi), (vo,Xi) and (y Q ,yi). Correspondingly 
we have four equations with the four unknown functions 
Po(w\x Q ), P {w\y ), Pi(tu|xi) and Pi(iZ>||/i). These equa- 
tions are the row elements of the vector equation 



RP = S, 



(A6) 



where 



R 



( r(Q\xQ,xx) r(l|x ,a;i) 

r(0|xo,j/i) r{l\x ,y\) 

r(0\y ,xi) r(l\y 0) xi) 
r(0\y o , Vl ) r(l\y , yi ), 

(A7) 



P 



and 



P {w\x ) 

Po(w\y ) 
Pi(w\ Xl ) 
Pi(w\ yi ) 



( 1 + w ■ v(xq, X\) 
l + w ■ v(xo,yi) 
1 + w ■ v(yo,xi) 
\ 1 + w ■ v(y ,yi) 



(A8) 



(A9) 



If R would be invertible, P n (w\x n ) would be a linear 
function of w, but this is impossible. It can be proved by 
contradiction. Suppose that the conditional probabilities 
are linear functions of w, that is, P n (w\x n ) = C n (x n ) + 
d n {xn) • w. They have to be non-negative and smaller 
than or equal to 1. It is clear that P n (w\x n ) are equal to 1 
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at most for only one vector w, but this is in contradiction 
with the fact that the conditional probabilities have to 
be equal to 1 in the one-dimensional manifold of vectors 
w such that f n (w) — x n . Thus, the equation 



det-R(x ,a;i, yo,Vi) = 



(A10) 



has to be satisfied for any x n and y n . Its solution is found 
by differentiating in j/q arL d y\, 



d 2 



dyodyi 
that is, 



det R(x ,xi,y ,yi) 



= 0, (All) 



yo=xo, yi=xi 



r(0\x ,xi) 2 r(l\x ,xi) 2 d Xo d Xl lo{ 



r(l|x ,xi) 
r(0\x ,xi) 



= 0. 



The functions r(n\x) cannot be identically equal to zero, 
thus the only acceptable solution that satisfies also the 
constraint (IA2D is 



r(0|a;o,xi) 
r(l\x ,xi) = 



fco(x ) + fci(xi) 

kp{xo) 
ko(xo) + ki(xi) ' 



(A13) 
(A14) 



where ki(Xi) are positive functions of the variable Xi. 
It is easy to check that this is also a solution of 
Eq. (|A10p . Thus, the fact that the conditional proba- 
bilities P n (w\x n ) cannot be linear in w allows us to find 
a constraint for the probabilities r(n\x). 

Since the determinant of R is equal to zero, there is a 
constraint also for the S because of Eq. (|A6|) . Let the row 
vector u T be the left eigenvector of R with eigenvalue 0, 
then we have from Eq. (|A6[) that 



(A15) 



•5 = 0. 



By means of Eqs. (|A13IA14p . after a bit of calculations 
we find that the left eigenvector of R with zero eigenvalue 
is 



M^o) -1 + fcitzi) -1 

-koixo)- 1 - /ci(yi) -1 
-koiyoy 1 - fciOi)- 1 
fcoM" 1 + fci(2/i)~ 1 



(A16) 



From Eqs. (|A9IA15IA16[) we obtain by differentiation the 
condition 



9 vo,vi (f T ■ S) 



■'ya.yi 



i , l 







v(x ,xi) 



= 0, 



(A17) 



that is satisfied if 



v(x ,Xx) 



k (x ) h(xi) 



[go(x ) +gi(xi)] , 

(A18) 



where go(%o) an d gi(%i) are generic vectorial functions. 
This is the inverse of the equation x — f(v). The func- 
tions g n (x n ) and k n {x n ) are constrained by the equation 
v 2 = 1, that is, 



[<7o (so) + gi(x 1 )} 2 = 



1 



k (x ) ki(xi) 



(A19) 



We will return to it later on. 

From Eqs. (|A5IA13IA14IA18|) and the identity S 
[1 + w ■ v(x)] /2 we find that 



E 



P n {w\x n ) - - 



1 



k^ixn) - 7;w-g 



0. (A20) 



Note that each term of the summation depends only on 
one of the variables xo and x\. Thus, the conditional 
probabilities have the form 



P (w\x ) = k (x ) 
Pi(w\xi) = fci(xi) 



~w ■ 9o{xo) + H(w) 
-iS-gi(xt) -H(w) 



\ (A22) 



where H(w) is an additional function independent of x. 
Since the probability Pq(w\xq) has to be equal to 1 for 
w = v , we find that 



H(v) 



1 



k (x ) 



v-go(x Q ) 



(A23) 



where Xq = fo(v)- Similarly, from the condition 
-Po^l^o) = 1 we obtain the equation 



H(v) = —- 



ki(xi) 



v-g\{x x ) 



(A24) 



This last equation can be derived by Eqs. (|A18IA23I) and 
the constraint v 2 = 1. It is interesting to note that the 
model we are constructing works only for a subset of 
preparation states v, as shown in Sec. IIII1 thus the func- 
tion H(w) is not necessarily given by Eqs. (|A23IA24I) if 
w is outside that subset. 

At this point we have almost everything, the last step 
is to find the functions k n (x n ) and g n (x n ) that solve 
Eq. (|A19|) . Differentiating this equation with respect to 
xq and x\, we have that 



dg {x ) dgi(xi 



d 



l 



( d 



l 



dxo dxi \dxo ko(xo) J \dxi k\(xi) 1 

' (A25) 

that is, the Minkowski inner product between the two 
four-vectors 



a(x ) 



go(x ) + io , , 
i i Yn ; p(si) 



Si(si) +7i 



(A26) 
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is equal to 0, 7„ and Xn being constant vectors and 
scalars, respectively. Using the Einstein notation, on the 
index contraction, the constraint is 



a fl (x )p»(x 1 )=0. 



(A27) 



It is important to note that a and /3 depend only on 
one of the variables xq and x\. Using Eq. (|A27|) . Equa- 
tion (|A19|) becomes (50-71) 2 - \h ~ Xi) -(7o + 7i) 2 + 

(xo+Xi) 2 = -7o) 2 + (± -Xo) 2 -The left-hand and 
right-hand sides depend only on xq and xx, respectively, 
thus they have to be equal to a constant rg, 



(.90 - 7i) 2 - 
(5i - 7o) 2 - (£ 



Xi 



Xo 



(70 +7i) 2 - (Xo +Xi) 2 
-To- 



rn 



(A28) 



These equations and Eq. (|A27I) are a convenient resettle- 
ment of Eq. (|A19|) . 

It is interesting to observe that the conditional proba- 
bilities P n (w\x n ) and Eqs. (|A18|A19[) are invariant with 
respect to the transformation 



9o 



9o + t,gi -> gi - t, 



(A29) 



for a generic vector t. This means that the constant 
vectors 70 and 71 are redundant and for example one 
could set 71 = 0. 

In order to solve Eq. (|A27|) . we have to consider two 
possibilities: (a) one of the two four-vectors a and /3 
spans a one-dimensional vectorial subspace and the other 
one lives in the orthogonal subspace; (b) the two vectors 
span orthogonal two-dimensional spaces. In the former 
case no solution exists such that v(xq,xi) is locally in- 
vertible (only a one-dimensional subspace of the Bloch 
sphere is represented). Thus, we consider the latter case. 
The most general solution up to rotations, transforma- 
tion (|A29|) and variable change x n —> F n (x n ) is 

COS Xq 

/7< 1 1 ' 0) = ( sin x sin j , = cos 9 cos x + s 

COS 9q CSC X\ 





cotxi sin#o 



CSCXi — S, 



(A30) 

where 9 and s are two free parameters. Only the former 
is present in the mapping x — > v(x), given by Eq. ([A18|) . 
We have that 



v(xq,Xx) 



1 



1 + cos 9qu x 



cos #0 + u x 
sin 9qu v 
sin 9qu z 



(A31) 



where 



u x — smii cosxo, 
u y = sinxi sin^o, 
u z = cosx\. 



(A32) 



The mapping is bijective by removing the poles at X\ 
and x\ — 7r, 



< x Q < 2tt, 

< Xl < IT. 



(A33) 




FIG. 2: Two orthogonal coordinate systems for 80 = 0.5 rad 
(left) and 80 — 1 rad (right). 



The inverse of the vectorial function v(x) in Eq. (|A31|) 
gives the functions f n (v) = x n , whose trigonometric 
functions are 



sinio = v y sin 6*0/ A, 
coscco = (v x — cos6'o)/A, 
smxi = coseo , 
cosxi = /- sinfl " , 

^ 1 — V x COS Wo 



(A34) 



Wy sin 2 1 



with A = (w x — cos ^o) 

Equation (|A31I) provides a set of orthogonal coordinate 
systems of the sphere. The spherical coordinate system 
is obtained with 9 — ir/2. Each system is mapped to 
another one by means of the Mobius transformation [3] ■ 
In Fig.[5]we report two coordinate systems for 9o = 0.5, 1. 
Both of them have two poles, but with different angular 
distance, that is equal to 29q. 

From Eqs. (|A13IA14IA30|) we have that 



\x) 



sin x\ (s-fcos #o cos xq) 

1+COS #o COS Xq sin Xi 

r(l\x) = . 

^ 1 ! l+COS Po COS Xq Sill X\ 



They are positive if 

cos#o| < s < 1. 
H(w) is obtained by Eq. (|Al3| . 

s- A 



H(w) 



(A35) 
(A36) 

(A37) 
(A38) 



Finally, the conditional probabilities for the events are 
given by Eqs. (IA21IA22|) . 



P (w\x) 



(w x — cos 9q) cos x + w y sin x sin 9q 



2(s + cos 6*o cos a;) 



(A39) 
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P (,TiM - i j- ( Wx C0S 6 ° Z jj + Wz cos x sin ^ + A sin x The model in Sec. [TTT] is obtained for 6> = tt/2 and 

2(1 — s smx) s=l. 

(A40) 
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